Self-energy of strongly interacting Fermions in Medium: a Holographic Approach 



Yunseok Seo^, Sang- Jin Sin''' and Yang Zhou^* 
t Department of Physics, Hanyang University, Seoul 133-791, Korea 
* Center for Quantum Spacetime, Sogang University, Seoul 121-742, Korea 

We consider the self-energy of strongly interacting fermions in the medium using gauge/gravity 
duality of D4/D8 system. We study the mass generation of the thermal and/or dense medium and 
the collective excitation called plasmino, by considering the spectral function of fermion and its 
dispersion relation. Our results are very different from those of the hard thermal loop method: for 
zero density, there is no thermal mass or plasmino in any phase. Plasmino in deconfined phase is 
not allowed in D4/D8 set up. In the confined phase, there is plasmino modes only for a window of 
density. 



I. INTRODUCTION 

The study of fermion self-energy in medium has a long 
history due to its fundamental importance in studying 
electronic as well as nuclear matter system. When the ex- 
citations are strongly interacting, perturbative field the- 
ory method cannot give a reliable result since the di- 
agrams should be truncated to the ladder or rainbow 
types, which can not be justified in strong coupling. Fur- 
thermore in the presence of chemical potential, the lattice 
technique is not much useful due to the sign problem. 
Therefore it is worthwhile to utilize the gauge/gravity 
duality for this tantalizing problem. The gauge/gravity 
duality was used to study the fuzzy fermi surface [1] and 
the non-fermi liquid nature [2H1] of the strongly inter- 
acting system, for recent developments, we refer to [SHE]- 

The weakly interacting field theory (QED or QCD ) 
results for the fermion self energy in medium can be sum- 
marized by the existence of the plasmino mode [9] and 
thermal mass generation of order gT. Plasmino is a col- 
lective mode whose dispersion curve has a minimum at 
finite momentum. However, recent study of the ther- 
mal field theory jTOl E], by solving Schwinger-Dyson 
equation numerically, showed that the thermal mass is 
reduced as the coupling grows. It raises an interesting 
questions what happens to the thermal mass and more 
generally to the hard and soft momentum scales and mag- 
netic mass scale, T, gT, g 2 T respectively, in the strong 
coupling limit. 

In this letter we study the dispersion relation and ther- 
mal mass generation using gauge/gravity duality and 
report a feature of plasmino in strongly coupled sys- 
tem. We consider D4/D8 model [H] and turn on a 
fermion field on the flavor brane world-volume with fi- 
nite quark/baryon number density 13 15j. In deconfined 
phase, the fermions represent the fundamental strings, 
connecting the horizon of black branes and the probe D8. 
In the zero 't Hooft coupling limit, it is a bi-fundamental 
representation. In the large 't Hooft coupling limit, D4- 
branes disappear and the color index of the fermions 
should disappear also. In confined phase, to describe 
a baryon we need to introduce Witten's baryon vertex 
j!6j . N c strings are connecting it to the probe D8/D8. 
The dynamics of connecting open strings define that of 



the baryon vertex as well as the D8. Since each string 
gives fermionic mode, total vibrational dynamics should 
be described by the composite operator coming from the 
product of all N c fermions. This is the baryon in con- 
fined phase which is fermion if N c is odd. We replace 
this composite operator by a massive fermion field. So 
our construction for baryon is rather phenomenological. 
We consider just one flavor brane case mostly. Previ- 
ously fermions on a probe brane in adjoint representa- 
tion, called mesino, were studied in [THUS], which are 
different from ours. 

By solving the Dirac equations in each phase, we ob- 
tain dispersion relations for the fermionic excitations in 
medium. Our results show that for zero density, there 
is no thermal mass generation and no plasmino in any 
phase, which is sharply different from weakly coupled 
field theory result. To get plasmino in deconfined phase 
we need to add quark mass as well as quark density. In 
DA/D8 set up, the quark mass is not allowed so plasmino 
is forbidden in deconfined phase. In the confined phase, 
there is a plasmino mode only for a certain window of 
density. 



II. PLASMINO IN FIELD THEORY 

We begin our discussion by giving a brief discussion of 
fermionic collective excitations in a plasma. The fermion 
propagator is written as 



G(p) 
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where £ = 7^5^ is self-energy. The gauge invariant re- 
sult is available in the hard (high temperature) thermal 
loop approximation (HTL) where fermion mass m can be 
ignored since it is small compared with Tor /i. The most 
noticeable effect of the medium is effective mass genera- 
tion. The HTL result of effective mass is given by [15] 

m'j — \g 2 CF (t 2 + ^2^j , where Cf — 1 for electron and 

Cf = 4/3 for quark. There are two branches of disper- 
sion curves lj = lo± (p) whose asymptotic forms are given 
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by 



= TV 



p « m f : uJ±(p) ~ m/ ± -p , 



(2) 

p >> m/ : w±(p) ~ p . (3) 

w + is the normal branch and w_ is the one describing 
plasmino that has been extensively investigated [H IT9l - 
I3T] . For a review we refer to [351 [33J . The presence of 
plasmino is important since it may enhance production 
rate of the light di-lepton [30] . 

In HTL approximation, fermions can be regarded as 
massless therefore two branches u>± can be characterized 
by chirality and helicity. The ratios of chirality and hc- 
licity for uj± are ±1 respectively. In our work, we do not 
neglect the fermion mass so we consider the presence of 
plasmino as the presence of minimum of dispersion curve 
at finite momentum. We will conclude that the plasmino 
exists if the dispersion relation has a minimum at finite 
momentum. 

Notice that plasmino in HTL approximation exists in 
high temperature whatever the density is. However, in 
the strongly coupled limit, we will show that plasmino 
disappears at zero density in deconfined phase. And it 
can survive only for a certain window of chemical po- 
tential in confined medium. Especially, at zero chemical 
potential, there is no plasmino. We find that the con- 
stant value 1/3 in Eq.([2| will be replaced by a function 
of density. See figure 7. 

III. HOLOGRAPHIC SET UP 



Let us now set up a holographic model to calculate 
the fermion self-energy, first in confined phase. We use 
Sakai-Sugimoto (SS) model [ISj where probe D8/D8 is 
embedded in black DA branes background. To introduce 
finite density, we turn on J7(l) gauge field on the probe 
brane. The sources of the !7(1) gauge field are end points 
of strings which are emanating from horizon in decon- 
fined phase and from baryon vertex in confined phase. 

Confined Phase The geometry of confined phase with 
Euclidean signature is given by 

/tt\ 3/2 

ds 2 = [ — ] (5^dx>*dx u + f(U)dxj) 

3/2 / dU 2 
CJ \f(U) 
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where both the time and the Kaluza-Klein direction are 
periodic: x° ~ x° + Sx , X4 ~ £4 + 6x4 and f(U) — 

(Ait\ 2 R 3 



1 Ukk 

1 C/3 



Ukk — (x)" fx 1 ■ Here following the orig- 
inal Sakai-Sugimoto model, we consider only trivial em- 
bedding of the flavor eight brane where D8 and D8 brane 
are located at the antipodal position in X4 direction. The 
dynamics of probe branes with U(X) gauge field on it can 
be governed by the DBI action as follow 



Sd8 — —Tg 
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(5) 



where AT = T^il^v^Sx R 5 / g s , and is the unit four 
sphere volume and V3 is the three Euclidean space vol- 
ume. Notice that in the case of antipodal configuration 
of D8/D8, the embedding becomes trivial i.e. x'^ir) = 0. 
For later convenience, we use dimensionless quantities 

r = U/R, r = Ukk/R , a = 2va'A /R. (6) 

From the equation of motion for U(l) gauge field, we get 
conserved dimensionless quantity ra i r ) _ _ 

D , where D is an integral constant representing the 
baryon density. 

Chemical potential can be defined as the value of ao (r) 
at the infinity, once proper IR boundary condition is im- 
posed at r = tq . In confined case, the baryon vertex 
which is DA brane wrapping on S 4 with Nq fundamental 
strings plays the role of source of the U(l) gauge filed. 
The energy of the source should be included in the total 
free energy of the system and contributes to the chemical 
potential also. 

Accordingly we set ao(ro) = m/q to include the baryon 
mass to in the chemical potential fi: 



/1 = m/q + (Iq = m/q + / a' dr 



(7) 



In confined phase, po is the contribution only from the 
gauge potential and to is identified as baryon mass, which 
can be calculated from DBI action of the baryon vertex. 
It is related to the 5 dimensional Lagrangian mass TO5 
as we will see later. Notice that such identification is 
not true in bottom up AdS model, where bulk fermion 
mass should be related to the conformal dimension of an 
operator. 

Deconfined Phase The deconfined geometry is given 
by double Wick rotating x 4 and time from the confined 
one. The DBI action for probe brane can be written as 
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The gauge field profile on eight brane is solved as be- 
fore, -. ra °^ — jy where D' is a constant. In 

^Jr- 3 (l-a' a (ry) 

deconfined phase the end points of fundamental strings 
play the role of sources of the U(l) gauge filed. Since the 
probe branes always fall into the black hole horizon, the 
fermion representing quark is massless in this phase. 

The IR boundary condition of ao field determined 
by the regularity condition at the horizon is given by 
ao(rfr) = , where rjj = Uh/R ■ This condition is also 
consistent with the argument in previous section. In de- 
confined case, mass of source is zero hence it does not 
contribute to the IR boundary condition for chemical po- 
tential. 

The schematic configuration of probe branes are drawn 
in Figure [T] 
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(a) 



(b) 



FIG. 1. (a) Embedding of D8/D8 brane in confined phase. 
D4 baryon vertices locate at r = ro. (b) Embedding of probe 
branes in deconfined phase. The source of U(l) gauge field is 
located on the black hole horizon. 



IV. 



FERMIONIC SPECTRUM 



Now we study the fermionic spectral function in the 
holographic dual system. We put a probe fermion field 
in the world volume of flavor D8 brane in the Sakai- 
Sugimoto model. We ignore the <S* 4 part and work in the 
effective 5 dimensional world volume following the origi- 
nal work of Sakai-Sugimoto model. The induced metric 
is written as ds 2 = g^dx^dx" + g rr dr 2 . We use the 
minimal action 



S = 



d 5 xy/=g~ i (i>T M D M il> - m 5 ^V) 



(9) 



where the covariant derivative is Dm = dM + \^abM^ ab — 
iqAjrf ■ Here M denotes the bulk spacetime index while 
a, b denote the tangent space index. After a factorizing 
tp = (-ggrryi/ie-ivt+ikix'ig ; thc Dirac equation for * 
can be give by 



^~/g~(T r -d r - m 5 ^)V + iKJ&S = , (10) 



where K M = (-v(r),ki) , v(r) = (w + qa )^ '-gu/ 'gtt ■ 
Following the procedure in [3] , we rewrite the Dirac equa- 
tion ( 10 ) in terms of two component spinors 



(d r +m 5 y/g^a 3 )$ a = v Vr/Sn(^M r ) + (- i r fco ' 1 ) $ c« , 

where a 1 are Pauli matrices and a = 1,2 denoting the up 
and down two spinor respectively. Further decomposing 
$! = (yi,2i) T ,$ 2 = (y2,z 2 ) T , we get equations of 
motion for y a and z a . When a = 2 we have 



(d r + m 5 ^/g^)y 2 (r) = y/g rr /gu{+v(r) + k)z 2 (r) (12) 
(d r - m 5v fg^)z 2 (r) = \/g rr /gu(-v(r) + k)y 2 (r) (13) 

By replacing k by —k, we obtain the equations of motion 
for yi and z\. We would like to define the following new 
variables G\{r) := yi{r)/z\(r) ,G 2 (r) := y 2 {r)/z 2 (r) . 
The retarded Green function in boundary field theory is 
obtained as 



Gf 2 = lim e 



-SmsRr 1 



Gl,2(r)| r =i/ e 



(14) 



where G\ and G 2 satisfy the following equations 



-d r G a + 2m 5y ^iG a 

(~l) a k + v(r) + ((-1) Q - X fc + v(r)) G 2 a . (15) 



Now the remain task is to solve ( 15 ) by imposing proper 
boundary conditions. 

Confined phase In the confined phase of Sakai- 
Sugimoto model, v(r) function appearing in the Dirac 
equation is given by v(r) = uj + qao(r) . where the U(l) 
flux solution on flavor brane is obtained from DBI action 
as 



ao(r) = /i + 



dr 



m-'D 2 

f 5 + D 2 



1/2 



(16) 



Notice that g rr = i? 2 r _3 / 2 /(r) _1 which diverges at ro. 
The boundary conditions for flow equation ( 15 ) can be 
found by requiring ( 15 ) regular at r = r . They are given 
by 



G a (r ) = 



—mR + \/m 2 R 2 + fc 2 — d> 2 

(-l) a k-UJ 



(17) 



m . And we define a 4 dimensional vac- 
^by/9ii\r=r ■ Notice that imposing the 



where uj = uj 
uum mass m 

boundary condition for retarded (advanced) Green func- 
tion corresponds to oj — > ui + it (u) u) — ie). 

Deconfined phase In the deconfmed phase, v(r) func- 
tion is given by v(r) = "+^j( r ) anc | the electric flux is 
given by 



a Q (r) = n + 



dr 



D 



12 



D' 2 



1/2 



(18) 



Due to the black hole horizon, we imposing the in-falling 
boundary condition at r = Th- The boundary condition 
at r = rn is obtained as 



Gi, 2 {r H ) = i 



(19) 



The fermion dispersion relation can be obtained by 
searching for poles of spectral function, which is the imag- 
inary part of retarded Green function. We solve (151 nu- 



merically with IR boundary conditions (17) in confined 
phase and ( 19 ) in deconfined phase. 



V. NUMERICAL RESULTS 

Now we discuss the results of fermionic spectral func- 
tion. First in the deconfined phase, the self-energy term 
gets some imaginary part due to the in-falling bound- 
ary condition. If the imaginary part is not large we can 
read off the dispersion relation from the positions of the 
maxima of spectral function. The 3D plot of spectral 
function with zero density is shown in Figure [2] Since 
the fermions here are deconfined quarks, we set m = 0. 
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(a) (b) 

FIG. 2. Spectral functions of at zero density: the dis- 
persion curve passes through the origin and thermal mass 
vanishes. Top view and side view of 3D plot of ImGf with 
T = 1, m = 0. True range of both ui and k is [—5, 5]. 



<3 



G? 






/ — - — 










— *w=0 line 










G? 






0.20 








0.15 








0.10 
0.05 
0.00 
-0.05 
-0. 1 () 






r--q.Ha 






/ m = 0.1 
/,,-''' rt, = 1.22 


Nil 0.2 0.3 0.4 0Va6 0.7 











0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



k 

FIG. 3. Dispersion relations in confined phase. The upper and 
lower branches describe the normal fermion ui+ and plasmino 
uj- respectively. Dotted line denotes light cone. 



Our main question here is whether thermal mass can be 
generated by finite temperature effect in strong coupling 
limit. In Figure[2]the dispersion curve passes through the 
origin and this feature is independent of temperature al- 
though it is illustrated for T = 1. As a result, no thermal 
mass is generated and there is no plasmino in dcconfincd 
phase with zero density. The absence of thermal mass 
is actually one of most drastic difference compared with 
the weakly coupled field theory result. Namely 

vtlt — —/=gT in weak coupling, 
V6 

m T = in strong coupling. (20) 

The result of vanishing thermal mass is actually consis- 
tent with a recent claim made in by numerical study 
of Dyson-Schiwinger equation in the strong coupling re- 
gion. 

If we turn on finite density in the system, density ef- 
fect can generate effective mass. However for massless 
fermion we have an exact Green function Gr — i at 
k = 0, so it is impossible to observe peak structure near 
k = 0. We can not observe plasmino mode in finite den- 
sity for the massless fermion. 

What happen if we added a finite bulk fermion mass 
for curiosity? We find that density effect can generate 
effective mass as well as plasmino mode for large enough 
chemical potential. 

Now we turn to the confined phase. In this phase, 
all the peaks of spectral function are delta function-like 
peaks. In general the Fermi momentum is defined by 
Lo{kp) = . We plot dispersion relations cj = Lu±(k) in 
Figure [3] 

Observation of plasmino We observe a plasmino dis- 
persion relation characterized by the presence of the min- 
imal energy at finite momentum. As we change chemical 
potential, the slope of cj_ at k = changes. We plot 
the slope a(no) as a function of density in figure [5] In 
HTL approximation, the value of slope at k = is — | 
independent of density or temperature. 
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FIG. 4. Dispersion relations of G2 in confined phase. As 
chemical potential increases, dispersion curve moves down. 
Parameter m = 0.1. 
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FIG. 5. no dependence of a — 4r at k = 0. The curve is 
plotted only in the density window /ii < fj,o < fi2 where there 
is plasmino. Parameter m = 0.. Notice that a — 4r at k = 
has a constant value -1/3 for the weakly coupled field theory. 



The high density behavior of the dispersion curve is 
complex and we will report it elsewhere. We restrict our- 
selves to the density range where traditional plasmino 
mode exists. We could observe plasmino only in a chem- 
ical potential window n± < fiQ < /i2- When m = 0.1, 
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the window is given by \i\ = 0.69 , [i% — 1.94. As to 
increases, this window gets wider. Inside the window, as 
density becomes larger, dispersion curve moves down and 
bends more and more as shown in Figure [4] This may 
be compared with field theory result in weak coupling, 
where effective mass and plasmino are generated for any 
density. 

We can read off baryon mass in medium defined by 
to* := lu(0). We will show to* is a monotonically in- 
creasing function of chemical potential. The mass shift 
in medium is defined as 5m := to* — to. The normalized 
mass shift quantity x(tM)) '■— Sm/m is plotted in figure[6] 
Notice that for massless case, there is no mass correction 
at all. 




FIG. 6. a): fio dependence of \- b) : 1^0 dependence of kp. 
Parameter m = 0.1. 



If we turn off both charge and bulk mass in confined 
phase, we have exact solution G2O") = \J fz^j , which is 
independent of radial direction. 

Fitting dispersion relation Now we try to fit the dis- 
persion curve of plasmino. We write the Green function 
as 



G*(u,k) 



W - Sm - k+B 



(21) 



Although the numerator Z is some undetermined func- 
tion, the pole information is assumed to be contained 
only in the denominator. Parameter B and C can be 
written in terms of Sm = to* — to, Uf and Vf as follows 



B= k 2 F {l-v F ) c = 
dm + kpVF ' 



+ 2Smkp + kpvp 
Sm + kpVF 



(22) 



where Sm is mass shift, kp is Fermi momentum and vp 
is Fermi velocity, vp is defined as the slope at k = kp. 



Expanding (21) near Fermi momentum, we get 
G*(w,*)~ Z 



vp(k — kp) — £ ' 



(23) 



where the self-energy near Fermi momentum is obtained 

as 

E = Sm+ (1 - VF )( Sm + *f Vf) .(k-k F r+Q((k-k F n 



kp{Sm 



The vanishing mass shift limit Sm = corresponds to 
the massless fermion in confined phase. In this limit, 
self-energy becomes 

S = (1 - ^ F)VF ■ (fc - k F f + 0{(k - kp) 3 ). (25) 



VI. 



CONCLUSION AND DISCUSSION 



In this paper we discuss our observations on the char- 
acters of fermion's self-energy in the dense medium. By 
using gauge/gravity dual, we recover normal and plas- 
mino branch. In deconfined phase, we showed that for 
zero density, there is no thermal mass generation. In the 
weakly coupled thermal gauge theory, T, gT and g 2 T pro- 
vide three well separated scales, with physical interpre- 
tations: hard momentum, thermal mass for fermion(or 
electric screening mass for gluon) , magnetic mass respec- 
tively. These masses play the role of the infrared regula- 
tor of different scales. Such separation does not happen 
for large coupling ~ O(l), and we believe that we can not 
define any such physical scale either for strong coupling. 
In the confined phase, plasmino excitations are present 
only for a window of the density. The group velocity of 
the plasmino mode at zero momentum turns out to be 
density dependent rather than a constant, —1/3, show- 
ing the deviation of the HTL approximation. It is also 
worthwhile to notice that without medium effect, there 
is no mass correction. We found a simple empirical for- 
mula for plasmino dispersion relations. We expect that 
the phenomena we discovered here is common in many 
holographic models, so that it is a universal character. 

There is an issue whether TO5 can be calculated in a way 
following [18] . where one starts with massless fermion 
in 10 dimension and performs the KK reduction from 
10 dim to 5 dim. The answer is negative because D4 
background does not have a direct product structure. To 
see this, let us consider 



{~f.D x + -f.Dy - m w )$(x, y) = , 



(26) 



(24) 



where x represents the non-comact directions and y does 
the compact directions. For D3 background, where man- 
ifold is direct product AdS^ x S 5 , the "f.D y term produces 
~ (I + 2)/ R ac i s which can be interpreted as a mass term. 
While in our model, D4 background is not a direct prod- 
uct. Therefore the KK reduction in our case generates 
term ~ (Z + 2)/r which is a potential rather than a mass. 
One may call such potential term as a "mass function" 
which depends on the radius. We will study the effect of 
mass function in future work. 

Note added: In the closing period of this work, an ob- 
servation on two branches of dispersion relation appeared 
in fM] with emphasis on the spin physics. 
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